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ABSTRACT 

We compare the peculiar velocity field within 65 h" 1 Mpc predicted from 2MASS photometry and 
public rcdshift data to three independent peculiar velocity surveys based on type la supernovae, surface 
brightness fluctuations in ellipticals, and Tully-Fisher distances to spirals. The three peculiar velocity 
samples are each in good agreement with the predicted velocities and produce consistent results for 
Pk = fi^ 6 /bj(- Taken together the best fit (3k — 0.49 ±0.04. We explore the effects of morphology on 
the determination of (3 by splitting the 2MASS sample into E+SO and S+Irr density fields and find 
both samples are equally good tracers of the underlying dark matter distribution, but that early-types 
are more clustered by a relative factor 6e/&s ~ 1-6. The density fluctuations of 2MASS galaxies in 
8/i _1 Mpc spheres in the local volume is found to be cts.k = 0.9. From this result and our value of 
Pk, we find <T8(O m /0.3) a6 = 0.91 ± 0.12. This is in excellent agreement with results from the IRAS 
rcdshift surveys, as well as other cosmological probes. Combining the 2MASS and IRAS peculiar 
velocity results yields cr 8 (O m /0.3) - 6 = 0.85 ± 0.05. 
Subject headings: 



1. INTRODUCTION 

Peculiar velocities are a unique probe of the distribu- 
tion of mass in the nearby universe. The velocity of an 
object, such as a galaxy, is the sum of two contribu- 
tions: the cosmological expansion and the peculiar ve- 
locity, which arises from gravitational attractions from 
surrounding overdensities, which are dominated by dark 
matter. 

In the linear regime, the peculiar velocity v(r) is given 
by 

v(r) = ^/^ ro(r0 ^). (1) 



where S m (r) = {p — p) /p, and p is the average density 
of the universe. Equation (JIJ is a byproduct of the as- 
sumption that structure forms as a result of the growth 
of small inhomogencitics in the initial density field. For- 
mally, it is valid only in the linear regime where 8 < 1 
i.e., on scales larger than ~ 5 Mpc. With an all-sky mass 
density field, the resulting velocity field can be predicted 
and compared to observations. Note that this velocity- 
velocity (v-v) comparison method requires a complete 
description of the density field, in order to use Equation 

There are two common approaches to estimate <5 m . 
One approach is to assume simple parametric infall mod- 
els (e.g. "Virgo infall"), and fit for parameters of the 
model. The second approach, adopted here, is to assume 
that galaxies are tracers of the mass density field. In this 
context, it has become a common practice to employ 
the simplifying assumption of linear biasing, in which 
S g = bS m , where S g is the galaxy density contrast and b 
is the bias factor relating the mass-tracer (galaxy) fluctu- 
ations with the mass fluctuation field. With the inclusion 
of linear biasing, there is a relationship between two mea- 
surable quantities, S g and v(r), in terms of one unknown, 
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(3 = fl^f /b. The velocity predictions can then be com- 
pared to measured velocities obtained through secondary 
distance indicators to determine the quantity (3. 

In practice, one can use Equation Q to determine f3 
by several methods. A direct comparison of the well- 
known Local Group (LG) velocity with the predicted LG 
velocity, v(0), is one possibility. However, in practice, 
for any redshift survey the integral in Equation is 
limited at some distance i? max - Thus determinations of 
(3 by this method depend on the assumptions about the 
convergence of the dipole at the survey limit. 

If we approximate the contributions from beyond i? m ax 
as a dipole U, then we may express Equation Q as 



V(r) = 47 



d*r% CO ^4 



U. 



(2) 



Clearly a degeneracy exists between (3 and U for the pe- 
culiar velocity of the LG. By using peculiar velocity data 
for many objects, one can fit for both (3 and U. Alter- 
natively, by making v-v comparisons in the LG frame, 
vlq (r) = v(r) — v(0), then the dipole contribution U 
cancels out. Essentially what is being measured is infall 
into overdensities with known S g . 

It should be noted however, measured velocities are 
radial which means that our predictions, although three 
dimensional, will be converted to radial velocities. We 
chose to describe radial velocities in the LG frame, and 
calculate them according to 



ulg (r) = [v(r) - v(0)] 



(3) 



In this paper, we have reconstructed the galaxy den- 
sity field using the Two Micron All Sky Survey (2MASS) 
rcdshift data. The sample is selected, weighted and 
smoothed in order to make predictions of the local pecu- 
liar velocity field. We employ the VELMOD technique 
(Willick et al. 1997) to perform a maximum likelihood 
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analysis which compares the redshift and distance esti- 
mates to the derived gravity field. In Section we give 
the details of the selection, correction and reconstruc- 
tion procedures that have been applied to the redshift 
data. Section ?? describes v-v comparisons for the SN 
la, SBF and SFI peculiar velocity datasets. In Section 
0]we explore how splitting the density field by morpho- 
logical type affects our results. We repeat the tests of 
Section ?? using the morphologically-segregated fields to 
make our velocity predictions. In Sections 5 and 6 we 
present a discussion of the results and show comparisons 
to values in the current literature. 

2. 2MASS GRAVITY FIELD 

We now turn our attention to the construction of grav- 
ity field from the 2MASS redshift survey. In what fol- 
lows, we discuss the selection of the sample, general fea- 
tures and completeness corrections applied to the data. 
We discuss the properties of the luminosity function and 
how it relates to the weighting scheme used to define the 
redshift space density field. Finally, we discuss the tech- 
nique used to transform the data from redshift to real 
space. 

2.1. 2MASS Data 

The 2MASS dataset provides an all-sky view of the 
nearby galaxy population in the J, H and K s bands. 
Near-infrared light has several advantages: first, it sam- 
ples the old stellar population, and hence the bulk of the 
stellar mass, and second it is minimally affected by dust 
in the Galactic Plane. Details of how extended sources 
are identified are given by Jarrett et al. (2000). We use 
2MASS data to define a K s magnitude limited sample 
of galaxies. Redshifts and morphological type designa- 
tions (de Vaucouleurs et al. 1991) are drawn from the 
HyperLeda 2 database. 

The selection criteria used to define the sample are as 
follows: 

1. K20 apparent magnitude < 10.5 

2. Galactic latitudes |6| > 11.5° 

3. Distance r < 6500 km s _1 . 

The K20 apparent magnitudes are those magnitudes 
defined to be within the same circular aperture for a K s - 
band isophote of 20 mag arcsec -2 . The K20 apparent 
magnitude limit of 10.5 is chosen to yield a high redshift 
completeness (~ 90%) for the sample. 

The cut in galactic latitude is chosen to reduce the in- 
completeness associated with galactic extinction and con- 
fusion due to the dust and stars, which increases closer 
to the plane of the Milky Way (6 = 0°). 

The final selection criterion limits the volume to 6500 
km s _1 . Given our JT s -band magnitude limit, beyond 
this distance only L* or brighter galaxies appear in the 
sample, hence shot noise becomes large. Furthermore, 
previous work (Rowan-Robinson et al. 2000; Kocevski 
et al. 2004) has shown that there are few significant 
attractors beyond 6000 km s _1 , at least until the Shapley 
Concentration is reached at 15000 km s _1 . 
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After making the distance and latitude cuts, the final 
magnitude limited 2MASS sample comprises 5032 galax- 
ies. For each object in the sample, all of the required 
information is known (except for 152 with unknown mor- 
phological type). We have also cross-referenced these 
galaxies in HyperLeda to obtain group identifications 
(Garcia et al. 1993). This was successful for 25% of 
our sample, leaving 1815 galaxies with group designa- 
tions. Note that HyperLeda grouping information only 
extends to a velocity distance of 5500 km s -1 and has 
not been revisited even though additional galaxies have 
been added to the HyperLeda database. 

2.1.1. 2MASS Redshift Completeness 

The parent 2MASS K w < 10.5, |6| > 11.5° catalog 
has a redshift completeness of 90%. The completeness 
is not uniform across the sky, however, and therefore 
small corrections are needed. Since the redshift data are 
compiled from public sources, there are no clear bound- 
aries. We note, however, that redshift completeness is 
likely affected by two selection biases: location of obser- 
vatories in the Northern or Southern Hemispheres, and 
extinction in the Zone of Avoidance (ZoA) at low Galac- 
tic latitudes. Note that while the sample here is K-band 
selected, and hence is much less affected by extinction, 
redshifts are typically obtained in the optical and so the 
redshift completeness will still be affected by extinction. 

Fig. shows the 2MASS redshift completeness per- 
centages in terms of the galactic latitude and declina- 
tion (D). While the redshifts are highly complete at high 
Galactic latitudes, there is a decline in the complete- 
ness at lower |6|. In the Northern celestial hemisphere, 
the completeness is very high, but drops in the southern 
hemisphere. We also find that the redshift incomplete- 
ness is not correlated to the magnitude. 

Since the overall completeness is high, and hence the 
corrections are small, we model the incompleteness with 
a simple functional form. We assume that the probability 
of observing a galaxy in our sample with a redshift is a 
separable function P cz (b,B)= P X (&)P 2 (D). 

Fig^shows a completeness is ~ 1 at high Galactic lati- 
tudes, with a decrease completeness towards the Galactic 
plane. After trial and error, for Pi(b) we adopt the func- 
tional form 



Pi(b> 11.5°) = l-exp 
Pi(6< -11.5°) = l-exp 



sin b 
sin 6n 
sin 6 



sin 6s 



(4) 



For declination, the northern hemisphere is quite com- 
plete whereas the southern hemisphere is less so. Fig, 
^ shows no strong dependence on declination within a 
given hemisphere, so we have chosen to simply fit with a 
step function 



P 2 (D>0)=D N = 1 
P 2 (D <0)=D S 



(5) 



The likelihood of the given set of galaxies with and 
without redshift is 



C cz = n p i(M p 2(D0 x J[ 1 - [Pi(6«)P 2 (D J )] (6) 
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Fig. 1. — Left: 2MASS redshift completeness is plotted as a function of galactic latitude. The high latitude regions are very complete, 
except for b ^ —90° values where the observed low value is due to a very small number of total galaxies. As b — > 0° there is a decrease 
in the fraction of measured redshifts due to the difficulty of obtaining redshifts close to the galactic plane. The area |6| < 10° is excluded. 
Right: Redshift completeness in equal area bins of declination for sin|fej > 0.2. Note the higher incompleteness in the southern equatorial 
hemisphere. 



where the first product sum is over galaxies with mea- 
sured redshift and the second is over those having been 
detected with no associated redshift. 

The maximum likelihood solution yields best fit pa- 
rameters 6n = 8.33° and 6s = —11.9°, and Dn = 1 
and Ds = 0.933. We correct for incompleteness in b and 
declination by weighting the galaxies with redshifts by 
l/[Px(6)P 2 (D)]. 

2.1.2. Cloning the Zone of Avoidance 

Due to increasing redshift incompleteness at low galac- 
tic latitudes, only galaxies satisfying |6| > 11.5° were 
selected. Leaving the unsurveyed regions empty would 
cause systematic errors in the dynamical model, because 
this region would behave as a void. This would have the 
effect of creating a false outflow, as can be understood 
through the 8 term in Equation |JT]|, 

There are several ways in which this effect can be cor- 
rected. For instance, the volume can be filled with a 
uniform distribution of galaxies that exactly match the 
surveys' average density. This method has the benefit 
of alleviating the systematic errors associated with the 
empty space, and since 6 = in this region, there is no 
gravitational effect. 

We chose instead to fill the plane by interpolating or 
cloning adjacent, equal area regions above and below the 
galactic equator (Hudson 1993). Due to the infrared 
nature of the survey, a galactic latitude cut as low as 
|&| = 11.5° could be made, which meant only a relatively 
small strip had to be cloned, as opposed to optical sur- 
veys which typically have higher limiting |6| (Santiago et 
al. 1995; Giuricin et al. 2000). Hudson (1994) found 
that if uniform density was assumed then the values of 
(3 are only ~ 10% higher. However, surveys in the ZoA 
(Staveley-Smith et al. 2000) suggest that cloning is a 
better approximation. 

In our cloning procedure, the K20 apparent magnitude 
and redshift, along with all of the other associated prop- 
erties of a cloned galaxy remain the same as the parent, 



with the exception of a new galactic latitude which de- 
pends on the parent and cutoff galactic latitudes through, 
sin b clonc = sin b parent - sin b min . 

2.2. Luminosity Function: $ (L) 

In order to correct for galaxy incompleteness due to the 
imposed flux limit, it is necessary to know the luminosity 
function (hereafter LF) . The empirical LF is often fitted 
by an analytical expression first described by Schcchtcr 
(1976), 





L ' 


dL 


(f) exp 


~L*_ 


~L* 



in which L* (or equivalently M* if we describe it in terms 
of absolute magnitudes) is the fiducial magnitude and 
represents the point in which the LF changes from a 
power law (with slope a) to an exponential. 

To estimate the LF parameters, we adopt the density- 
independent method of Sandage et al. (1979). In Sec- 
tion ^. 3. 31 below, we will outline the redshift transforma- 
tion method which uses the smoothed, weighted density 
field to make iterative corrections to the distance of each 
galaxy based on an assumed (3. After each iteration, 
when the distances to every galaxy have been updated, 
we recalculate the LF parameters. Thus in contrast to, 
for example, deep samples in which redshift is used as a 
proxy for distance, in our case a self-consistent set of LF 
parameters for each value of /3 is derived. Nevertheless, 
at low redshifts the uncertainty in a galaxy's distance 
due to its peculiar velocity is large, so we have used only 
galaxies in the distance range 1500 km s _1 to 6500 km 
s" 1 to estimate the LF parameters, and excluded faint 
galaxies with Mk 20 > —18. 

We obtain best fitting values Mt = -23.42, a — 

-0.982 and <j)* Kw = 0.011 (in h 3 Mpc~ 3 ) from the analysis 
of the real space galaxy distribution (at the best fitting 
P, see Section I3.3|l . This is in agreement with the re- 
sults of Cole et al. (2001) who used the 2MASS and 2dF 
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galaxy redshift surveys to derive values for the parame- 
ters — 51og/i , a and $k in the K s band yielding, 
-23.44±0.03, -0.96±0.05 and 0.0108±0.0016/i 3 Mpc" 1 
respectively. The likelihood analysis of the LF was per- 
formed with and without the inclusion of cloned galaxies. 
We find that there is no appreciable difference in either 
case, as expected since the Sandage et al. method is 
density-independent. In Section 0] we explore how mor- 
phology affects the LF parameters by segregating the 
density field into early and late- types. 

2.3. Transforming Redshifts to Real Space Distances 

2.3.1. Weighting the Sample 

Knowing the LF and completeness of the sample, we 
can calculate the selection function (the probability that 
a galaxy is included in the redshift sample) at any dis- 
tance r as, 



<t>(r)i 



J4-jrr 2 /„- 



$ (L) dL 



1 



r > r n 



r < r n 



(8) 



where / m ; n is the flux corresponding to -ft^ojim = 10.5. 
In practice the luminosity function is poorly constrained 
at the faint end, and hence we place a lower limit on 
the integral at L min = 47rr^ in / min . L min corresponds 
to M K = K 2 q.\uu - 51og(r min /10pc), which is -18 for 
r m in = 5/i _1 Mpc. 

We assign extra weight, l/(f>(r), to galaxies with mea- 
sured redshifts to account for the galaxies in the same 
volume that were not detected or failed to meet the se- 
lection criteria. This is equivalent to adding galaxies at 
every position where we have a galaxy with a measured 
redshift. 

Weights are assigned to all galaxies in our sample limit, 
t < -Rmax- Galaxies outside of i? ma x are set to a weight 
of zero. The sum of the weights divided by the volume 
within i? max yields the average number density. 

2.3.2. Smoothing the Sample 

The galaxy distribution obtained from redshift surveys 
is a point process. If we were to apply linear theory 
directly, the predictions would diverge near the weighted 
points. We smooth the density field so that the velocity 
field is continuous and so that linear theory applies. This 
is accomplished by applying the following equation to the 
density field, 



«5(r) = 



P(*)-P = 1 
p n 



i 



w 



- 1 



( 9 ) 

in which W is a smoothing kernel and n is the number 
density of the N sa \ weighted objects. We have chosen 
to use a top-hat kernel in which W — |r — ^i\ 3 /r^ m for 
|r— Tj| < r sm and set to unity otherwise. We use a 
fixed smoothing length r sm of 500 km s _1 independent 
of distance. 

2.3.3. Reconstruction Procedure 

We now discuss the reconstruction method used to 
transform the measured redshifts into real space dis- 
tances, r, for the model. Given a smoothed all sky den- 



sity field, the peculiar velocity of any galaxy can be esti- 
mated via Equation J5J), and used to make distance cor- 
rections to the sample. We follow a similar method to 
that of Yahil et al. (1991) who use an iterative technique 
in which gravity is adiabatically "turned on" by increas- 
ing p. The following outlines a recipe for a self-consistent 
solution to real space density and velocity maps given 
only positions and redshifts. 

Prior to the start of the iterative procedure, we collapse 
the so-called "Fingers of God" (nonlinear redshifts dis- 
tortions resulting in an apparent radial stretching along 
the line of sight) by placing all grouped galaxies at a 
common median redshift, czlg> an d angular position. 

Initially, galaxies are assigned distances that are equal 
to their Local-Group-frame redshifts. A value of (3 is de- 
fined at each step of the iteration and the resulting pe- 
culiar velocity field is calculated. We perform 100 steps, 
defining a density and velocity field at each iteration from 
(3 = 0.01 to f3 = 1.00 in steps of 0.01. At each iteration 
the following steps are performed. 

1. The likelihood function for the LF is minimized to 
determine the best fit parameters and a number 
weight is assigned to each gala xy acc ording to its 
distance as described in Section \'Z. 3. II 

2. Galaxies are cloned to fill the void as described 
in Section f2 . 1 . 21 assigning each cloned galaxy the 
properties of its parent, including its weight. 

3. The average density n of all galaxies within i? max is 
calculated. This then defines the density contrasts 
S g . 

4. A peculiar velocity is calculated for each galaxy 
using Equation @. We average the current pecu- 
liar velocity with that of the past 5 iterations. As 
we step through /3, we are essentially turning on 
the gravity, and since increases in /3 are very small 
at each step, changes in weight and positions are 
also small over this range. The averaging has the 
benefit of damping out any unphysical oscillatory 
behavior 

With this information we update the distance of 
each galaxy according to the equation for total re- 
cession velocity, j-j = czLG,i — u(ri) (where i denotes 
the i th of iVg a i galaxies). Note that the distance to 
each galaxy in the sample changes at each itera- 
tion. In particular, after an iteration, a galaxy's 
distance may "move" it across R max . Therefore, at 
each iteration, we update the positions of all galax- 
ies with czlg < 9000 km s _1 , which is well beyond 
-Rmax- A galaxy residing in the region r > i? m ax 
and czlg < 9000 km s _1 is assigned a weight of 
zero, and therefore contributes nothing to the den- 
sity and predicted peculiar velocity fields. If a sub- 
sequent iteration brings it within i? ma x, it is as- 
signed a weight. Similarly, if an iteration brings a 
galaxy from r < i? max to r > i? m ax it is assigned a 
weight of zero. In this way, we maintain a density 
field based on a self-consistent distance limit -R m ax- 

We repeat the above steps at each iteration. Thus after 
the iterative procedure is complete, we have 100 density 
and velocity fields (specified by (3) 1 each of which can be 
compared to measured peculiar velocity data. 
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Fig. 2. — The isodensity contours for a slice through the Su- 
pergalactic Plane of the reconstructed 2MASS density field. The 
density field has been smoothed with a Gaussian kernel of 500 km 
s~ 1 . The heavy contour indicates mean density, 5 g = 0, and in- 
creases in steps of A5 g = 0.5. It can be seen that the zone of 
avoidance has been filled with cloned galaxies. Superimposed is 
the predicted peculiar velocity field for Pk = 0.49. 



2.3.4. Cosmography 

Fig. |21 shows the Supergalactic Plane maps of the re- 
constructed density field and the resulting velocity field 
for our best fit (3k = 0.49 as derived using the 2MASS 
rcdshift data. The galaxy density field is smoothed us- 
ing a Gaussian kernel with a smoothing radius of 500 km 

The plot shows the major structures that are located in 
the Supergalactic SGX-SGY Plane as well as the result- 
ing flow fields. Directly above the center (LG) in the pos- 
itive SGY direction is the Virgo Supercluster, with an as- 
sociated peak overdensity of S g = 2.7. At SGX — —4000 
to -3000 km s" 1 and SGY = 1000 to 2500 km s" 1 is the 
Great Attractor (GA). The peak in seen in Fig. [5]is coin- 
cident with the cluster Abell 3574 (r - 5300; 5 g = 3.4), 
but the highest peak in this region lies slightly off the 
Supergalactic Plane, and is coincident with the Centau- 
rus cluster at (r ~ 3600 km s _1 ; 5 g = 4.7) with a sec- 
ondary peak at Pavo II (Abell S0805; r ~ 4500 km s" 1 ; 
S g = 4.6). Note that what may be the most massive clus- 
ter in the GA region, Abell 3627 (Kraan-Korteweg et al. 
1996) is not in our sample because of its low Galactic 
latitude (b = —7°). However, the dynamical influence of 
a single cluster, even a massive one such as Abell 3627, is 
small. Instead the peculiar velocity field is more strongly 
affected by the large-scale overdensity of the supercluster 
as a whole. It is these large-scale overdensities that are 
crudely reproduced by our cloning procedure. 

The Perseus-Pisces (PP) supercluster shows the largest 
overdensity of 5 g = 5.5, coincident with the Perseus clus- 
ter, located near the sample limit and at the edge of 
the ZoA, in the SGX, -SGY direction. Note however 
that the overdensity may be enhanced due to the cloning 
of Perseus. The region at SGX = 2200 km s" 1 and 
SGY = —5000 km s _1 is harder to classify. It has an 



overdensity of S g = 2.6, and appears to coincide with 
Abell 168A/Abell 1194 overdensities. The Coma cluster 
lies outside of our sampled distance limit at around a red- 
shift of cz = 7000 km s -1 . In contrast, the -SGX, -SGY 
quadrant is almost completely devoid of overdensities, 
and is dominated by the Sculptor void which exhibits a 
weak velocity outflow. 

3. VELOCITY- VELOCITY COMPARISONS 

With the full velocity field modelled, we now make 
comparisons to peculiar velocity data sets. Below is a 
discussion of the properties of the peculiar velocity sets 
used, as well as the methods used in their comparison to 
the model. 

3.1. Peculiar Velocity Data 

We compare our predictions to three published peculiar 
velocity data sets. Each of which vary in sample size 
and typical distance errors and method of obtaining the 
secondary distance information. 

Our subsample of the Spiral Field I-Band (hereafter 
SFI; Haynes et al 1999a,b) survey consists of 836 galaxies 
of morphological type Sbc-Sc, in which distances have 
been derived from the I-Band Tully-Fisher relation over 
the full sky. Typical distance errors are on the order 
of 20% for our subsample, which extends to 5000 km 
s _1 . The characteristic depth of the sample, which is the 
weighted distance from which which the majority of the 
signal arises is 

4 = ^ (10) 

where the weights Wi = l/o~d,i and o~d,i is the distance 
error. For the SFI set, d c = 2700 km s _1 , so the charac- 
teristic distance error is o~d c = 540 km s _1 . 

The I-band Surface Brightness Fluctuation (SBF, 
Tonry et al. 2001) survey consists of 266 galaxies ex- 
tending to a distance of 4000 km s _1 . The uncertainties 
in the distance estimates varies inversely with the resolu- 
tion of the images. Typical distance errors for our sample 
are of the order of 8%, approximately half that of SFI. 
The characteristic depth is 1200 km s" 1 having ay c = 96 
km s _1 . 

The final data set is a compilation of supernovae of 
type la (SN la, Tonry et al. 2003) and consists of 59 SNe 
extending to a distance of 6000 km s _1 , with typical 
errors of 8%. The characteristic depth of this sample is 
2200 km s^ 1 with a Tc = 176 km s _1 . 

3.2. Method 

We apply two v-v comparison methods. Our primary 
method is a slight variant of the VELMOD method. For 
comparison, we also apply a simple and complementary 
X 2 fitting method. 

3.2.1. VELMOD 

The procedure outlined here differs slightly from the 
VELMOD method by Willick et al. (98). Our procedure 
uses two free parameters: /?, which scales the predicted 
peculiar velocities, and Hq, which allows for a rescaling 
of the published distances. In contrast, Willick et al. fit 
(3 as well as the parameters of the TF relation. 
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Specifically, we maximize the probability of each 
galaxy having its observed velocity, P(cz bs)- We con- 
struct the joint probability distribution of redshift and 
the (unobservable) true distance, 



P{cz ohs ,r) = P(cz ohs \r)P{r) 



(11) 



where the first term 



P(cz ohs \r) 



27TCT, 



■ exp 



[cz obs - (H r + /3u (r))Y 

2(7.2 



(12) 

is a description of the redshift-distance relation in terms 
of a Gaussian probability distribution. The predicted 
capped has three parts; an expansion velocity; a part due 
to the linear perturbations of the surrounding overdensi- 
ties, u(r), and a part which is often referred to as velocity 
noise, a v , assumed to be Gaussian, arising from strongly 
non-linear processes. For this calculation, o~ v is fixed 
at 200 km s _1 for field galaxies. This allows for uncer- 
tainties in the linear predictions, nonlinear motions and 
observational errors in cz Q ^ s . This value of a v is chosen 
becau se it gives reasonable reduced x 2 fits (see Section 
13.2.2(1 . although our recovered values of (3 are not sensi- 
tive to precise value chosen for o~ v . An additional velocity 
error was added in quadrature to a v to account for extra 
velocity dispersion for galaxies associated with clusters: 
o-ci = cto/a/1 + (r/r ) 2 , where ct = 700(400) km s" 1 
and ro = 2(1) Mpc are adopted for Virgo (Fornax) re- 
spectively (Blakeslee et al. 1999). 

The second term in eq. ITTl is the a priori probability 
of observing an object at distance r 



P (r) oc exp 



2a 2 . 



(i + <y 



(13) 



and is a Gaussian with mean equal to the estimated sec- 
ondary distance d with an additional term (1 + <5 g ), pro- 
portional to the number density of galaxies, to account 
for inhomogeneous Malmquist bias. The (1 + 5 g ) term 
is evaluated using the 2MASS density field, smoothed 
with a Gaussian filter of 500 km s _1 . In principle, for 
each peculiar velocity dataset, one should use <5 g for the 
appropriate morphological type, for example, 5e for the 
SBF dataset. In practice, as we will discuss in Section 
??, the correction arising from the (1 + <5 g ) term is small 
except for the SFI case. In Section 21 we show that the 
bias b for Spiral and later types appears to be very sim- 
ilar to that of the sample as a whole, thus justifying the 
simpler approach adopted here. 

The expression P(cz b s ) can be obtained by integrating 
over all possible distances, 



/>oo 

P{cz ohs ) = / P {cz ohs \r) P (r) dr . 
Jo 



(14) 



Equation (|14|) is then normalized over all possible veloc- 
ities. 

The VELMOD method then maximizes the log like- 
lihood, J2i \ n [P( cz i)], over all galaxies in the peculiar 
velocity data set as a function of the free parameters (3 
and Hq. 

3.2.2. x 2 method 



TABLE 1 

Summary of J3 determinations: VELMOD and x 2 

MINIMIZATION RESULTS. 



Sot: 


/^VELMOD 


13 2 


D.O.F 


x 2 


SN la 0.47 ±0.06 


u - 4o=t -0.05 


59 


42 


SBF 


0.46 ± 0.06 


n 44++0.05 
u - 44=t -0.07 


266 


238 


SFI 


0.55 ±0.05 




831 


756 


All 


0.49 ±0.04 









In addition to VELMOD, we perform a simple x 2 min- 
imization. The x 2 function is constructed from the dif- 
ferences between the observed cz b s and that which is 
predicted from a combination of our radial velocity and 
the galaxies measured distance. We minimize the func- 
tion, 



I C^obs {H d + /3u(d))Y 



(15) 



which is similar to that contained inside the parenthesis 
of Equation i|12fl . In this case however we do not evaluate 
it(r) at all possible values of r and marginalize, but in- 
stead use the peculiar velocity at the estimated distance 
d. To the extent that the peculiar velocity predictions 
do not change rapidly on the scale of the distance er- 
ror, a r , this approximation will be accurate. The error 
a cz — yJo- v 2 + ad 2 , where a v is the same as in Equa- 
tion 1(12(1 . Note that this application of the x 2 method 
neglects inhomogeneous Malmquist bias. Although it is 
possible to correct for the estimated distances for this 
bias (Hudson 1994a), we have not done so in this paper. 

3.3. Results 

The VELMOD analysis applied to the SN la, SBF 
and SFI peculiar velocity sets yields the likelihood con- 
fidence ellipses shown in Fig. U3 Table ^ summarizes 
the results for the peculiar velocity comparisons with 
those predicted by the reconstructed 2MASS density 
field. /Svelmod indicates the value obtained from the 
first moment of the marginalized (3 probability distribu- 
tion function (PDF) obtained from the VELMOD analy- 
sis. The /Jvelmod values are then combined to determine 
the best fitting (3 — 0.49 ±0.04, which is indicated as All 
in the last row of Tabled The (3 X 2 values are quoted for 
the x 2 minimization. There was only one free parameter 
in the x 2 fits corresponding to /?, and Hq was held con- 
stant at the best fit VELMOD value. Since there is only 
one free parameter in the fits, the number of degrees of 
freedom (DOF) is iV ga i - 1. 

With the adopted value of a v , the x 2 values are rea- 
sonable. The agreement between the SN la and the SBF 
samples is evident, and these two samples agree (within 
errors) with the slightly higher (3 determination for the 
SFI sample. All peculiar velocity datasets agree within 
the determined errors. The agreement is reinforced in 
Fig. 0] where we plot the measured SN la, SBF and SFI 
peculiar velocities against those predicted by 2MASS for 
each respective /?velmod ■ 

4. EFFECT OF MORPHOLOGY 
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0.2 



Beta 



Fig. 3. — The 68% and 90% joint probability contours on /3 and 
Ho are plotted for the VELMOD likelihood analysis for each pe- 
culiar velocity data set (SN la, SBF, SFI) and the derived 2MASS 
velocity field. Note that Hq is very close to unity in all cases. 
This joint distribution is then marginalized (over Hq) to obtain a 
probability distribution function for /3. 



It is likely that mass is related to light in a way more 
complicated than suggested by linear biasing. While such 
deviations may be small on linear scales probed by pe- 
culiar velocities, there remains the possibility of measur- 
ing such deviations. In particular, it is well-known that 
early-type galaxies are preferentially found in dense en- 
vironments. In the case of a toy model in which, for 
example, all mass of the Universe was in clusters, with 
no mass in the field, then one would expect the elliptical 
density field to be a better predictor of peculiar veloci- 
ties than that of spirals. So it is of interest to explore 
the predictions of density fields pre-selected by morpho- 
logical type. 

We have morphological types for nearly 97% of our 
sample. We consider two subsamples, an early-type 
(E+SO) sample having T < 0, and a late-type (S+Irr) 
sample having T > 0. Here T represents a numerical 
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Fig. 4. — Comparisons of SN la, SBF and SFI peculiar veloci- 
ties to those which have been derived from the 2MASS modelled 
velocity field out to a distance of 6500 km s —1 . Errors are defined 
by <?cz, but for SFI and SBF data only one in three error bars are 
plotted. Predictions are scaled to j3 = 1, so the slope of best fit 
corresponds to the value of /3. The solid line indicates the best fit 
to the individual data set shown, and the dotted line is the global 
fit Q3 = 0.49). 



code for the revised (de Vaucouleurs) morphological type. 
Our cloned sample contains -/Ve+so = 2486 galaxies. De- 
spite the -ftf s -band selection, the 2MASS rcdshift survey 
contains more late- type galaxies than early- types, such 
that iV s = 3795. 

We analyzed both in the same manner as outlined in 
Section 12.3.31 The LF parameters for the early-type 
galaxies are = —23.57 and a = —0.72 and the cor- 

responding values for the late-types are — —23.08 

and a — —0.79. Note that both of these are natter than 
the full density field, which may be a result of incom- 
pleteness in morphological types at the faint end. This 
trend is also observed in the recent determination of early 
and late-type LF's by Kochanek et al. (2001). They 



find M 
M 



K- 2 



= -23.53 ± 0.06 and a = -0.92 ± 0.10 and 
-22.98 ± 0.06 and a = -0.87 ± 0.09 for the 
early and late-type LF parameters respectively. They 
used 7 < K20 < 11.25 mag, with cz > 2000km s _1 and a 
limiting latitude of |6| > 20° for an approximately equal 
sample size. 

The derived density and velocity fields from the two 
morphological subsamples are shown in Fig. 03 The most 
distinct difference between the fields is that the early- 
types clearly reside in higher density regions. The density 
contrast is more significant in every clustered region, as 
is indicated by the increase in contour line density. Since 
the amount of clustering is quantified through the biasing 




Fig. 5. — The reconstructed E+SO (left) and S+Irr (right) density and velocity fields. Contours are as in Fig[5] Predicted peculiar 
velocities are for the best fit values from Table 2. While general features are consistent in both fields, the relative contrast indicates that 
the early-type galaxies are more clustered. 



TABLE 2 
j3 DEPENDENCY ON MORPHOLOGY. 



Set: 


/^VELMOD 




D.O.F 


x 2 


Late- Type (Spirals) 


SN la 0.46±+ u0 ^ 
SBF 0.45 ±0.06 
SFI 0.50±±g;|g 
All 0.48 ± 0.04 


o.48±ir b 7 

0.38±T°;°oI 
0.48±I°:°1 


58 
265 
830 


46 
244 
770 


Early-Type (E+S0) 


SN la 


0.32 ± 0.05 




58 


40 


SBF 


0.28 ±0.04 




265 


211 


SFI 


0.30 ±0.04 


n oo,+0.02 
u - oo=c -0.03 


830 


781 


All 


0.30 ±0.03 









parameter 6, we can also measure b via the (3 of best fit 
to the peculiar velocities. 

The results of the VELMOD and \ 2 minimizations for 
the morphological subsamples are given in Table |5J For 
the VELMOD determinations there is excellent agree- 
ment between the values of (3 derived from each of the 
peculiar velocity data sets for both early- and late-type 
density fields. There is also good agreement between 
the x 2 an d VELMOD techniques. The early- type den- 
sity field yields (3 values (~ 0.3) that are consistently 
lower than in the late-type case (~ 0.48). This suggests 
that the early-type density field is more strongly clus- 
tered (3s /Pe = bs/bs ~ 1.6 in the JsT s -band. 

In principle, it is possible to test whether early- or 
late-types are better tracers of the mass by comparing 
the quality of the fit from the two different density fields 
to the same peculiar velocity sample. For the SN la and 
SBF data sets, the peculiar velocity field predicted by 
early-types has a lower x 2 value, although this difference 
is only significant for the SN la sample. For the SFI data 
set, the opposite is true: the late-type density field yields 



a better goodness-of-fit. One might be concerned that, 
because the early-type density field is sparser (and hence 
more subject to shot noise), its predicted peculiar veloc- 
ities are noisier and hence the \ 2 mav be biased high. 
However, the early- type galaxies are also more strongly 
clustered, with the result is that the signal-to-noise ra- 
tios of the density fields are similar. To put it another 
way, the higher degree of noise in the density field is 
suppressed by the lower value of (3. We have performed 
tests to simulate the effects of sparse sampling on the 
goodness-of-fit statistic x 2 > ancl nn d no conclusive evi- 
dence for a difference between early and late- type density 
fields. However, further tests are required to understand 
the systematics in greater detail. 

5. DISCUSSION 
5.1. Potential Systematic Effects 

We now discuss the potential systematics which might 
affect our results. These include smoothing, external 
gravity contributions, sparseness of the sample, as well 
as the possibility of inhomogeneous Malmquist bias. 

Berlind at al. (2001) explored the effects of apply- 
ing various biasing schemes and smoothing lengths in 
the determination of (3. For the case of linear biasing, 
they find that the optimal smoothing length is roughly a 
3ft. _1 Mpc Gaussian sphere, which is close to our adopted 
5 ft^Mpc top-hat smoothing. In contrast, for example, 
had we adopted a top-hat smoothing of 8/i _1 Mpc, the 
results of Berlind et al. suggest that our derived value of 
f3 would have been biased too high by ~ 10%. Therefore 
our choice of smoothing length is close to optimal and 
we do not expect a significant bias in our determination 
of 13. 

The external velocity field can be reasonably described 
by the two main components (a dipole and a quadrupole, 
i.e. bulk and shear). Since the analysis was carried out 
in the LG frame, we need not worry about any external 
dipole contributions, however, this is not true for any 
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higher order contributions, such as the quadrupole, that 
may exist. Willick & Strauss (1998) suggested the exis- 
tence of a residual quadrupole velocity field that was not 
well modelled by the IRAS 1.2Jy gravity field. We find 
no evidence of a residual quadrupole in the SN la and 
SFI data. The SBF data suggests a marginally signifi- 
cant quadrupole, but the inclusion of such a term does 
not significantly affect the value of (3. 

Another possible systematic effect is related to the dis- 
creteness of the density fields used to make our velocity 
predictions. Since the galaxy density field is a sparsely 
sampled, this can lead to shot noise in the velocity pre- 
dictions. To test whether the sample size affects the de- 
termination of /3, we generated a sparse-sampled density 
field comprising of 50% of the full 2MASS density field. 
For this test, we find a small change (~ 2 — 3%) in the 
predicted value of (3, which is negligible compared to the 
random errors on /3. 

Finally, let us discuss the corrections made for inho- 
mogeneous Malmquist bias (1MB), which is a Malmquist 
bias arising from distance errors scattering peculiar ve- 
locity data away from overdensities. Since their redshifts 
are not scattered significantly, the resulting pattern mim- 
ics infall, leading to a biased estimate of (3. The strength 
of the bias increases as the square of the typical peculiar 
velocity error, o~d c , thus we expect the SFI sample to be 
most affected by the 1MB. The VELMOD technique in- 
cludes a (1 + S g ) correction to the distance probability 
function P(r) which corrects for this effect. It is never- 
theless of interest to estimate the effect of the 1MB on 
different samples, which can be done by repeating the 
analysis omitting the (1 + S g ) term. For the SN la and 
SBF samples, the correction is negligibly small (0.02) in 
j3. The SFI data set is more sensitive to 1MB; our result 
for this sample (3 = 0.55 (with 1MB correction) would 
have been (3 — 0.71 without the 1MB correction. For 
deeper samples with large errors, it is clearly important 
to account for 1MB. 

In summary, we have not identified any systematics 
which affect the result at a level greater than the random 
errors. 

5.2. Comparison with Other Results 

This paper is the first comparison between the 2MASS 
density field and peculiar velocity surveys. However, our 
results are consistent with previous analyses comparing 
the 2MASS dipole to the motion of the LG. Mailer et al. 
(2003) assumed n m = 0.27 and found b K = 1.06 ± 0.17 
which is equivalent to (3k — 0.48±0.08. Our result is also 
consistent with the upper limit (3k < 0.55 ± 0.2 derived 
by Erdogdu et al. (2005). 

TableOllists our results and recent (3 measurements ob- 
tained via velocity-velocity methods obtained using the 
IRAS 1.2 Jy and PSCz surveys. The 2MASS f3 values 
agree quite well with IRAS values. However, when com- 
paring the 2MASS predictions of (3 to the IRAS values, 
one should keep in mind that the IRAS samples are dom- 
inated by late-type galaxies. Consequently, the strength 
of the clustering, or bias, is likely to differ and the com- 
parisons of (3 may not be straight forward. With this 
in mind, it is preferable to transform the results into a 
form that is independent of the details of the density 
field used. This can be accomplished by noting that the 
linear biasing factor b relates the r.m.s. galaxy fluctua- 



tions, trg g, to the amplitude of r.m.s. mass fluctuations, 
erg, where the subscript 8 indicates that the density fluc- 
tuations have been averaged in top-hat spheres of radius 
8 /i _1 Mpc. If er 8ig is known, we may express our results 

as /3gCT 8 ,g = 

From the 2MASS K s band density field reconstructed 
in this paper, we have measured directly the r.m.s fluc- 
tuations in 8 Mpc spheres. Allowing for shot noise, 
we find a value of os,k — 0.90 (in reconstructed "real" 
space). We have not measured the uncertainty in this 
quantity, but note that this result is consistent with 
cr 8 K = 1.0 ± 0.1 found by Mailer et al. (2005) for 
a larger K s band 2MASS sample, although it is lower 
than the value found by Frith et al. (2005): 1.16 ± 0.04. 
Luminosity-dependent biasing may account for some of 
the difference between these results. Here we adopt 
a 8>K = 0.9 ± 0.1 to obtain the result cr 8 (f2 m /0.3) a6 = 
0.91 ± 0.12. We note that the is dominant contribution 
to the uncertainty is from the uncertainty in o~s,k- We 
expect that ct 8 k will be better determined in the near 
future from the 2MASS Redshift Survey (Huchra 2000). 

For the IRAS survey, Hamilton & Tegmark (2002) 
found that a 8 ,i = 0.80 ± 0.05. We quote ct 8 (an/0.3) - 6 
values in the third column of Table For the IRAS 
case, the v-v results have been combined to yield an av- 
erage IRAS value cr 8 (fi m /0.3) - 6 = 0.83 ± 0.06. There 
is good agreement between the various methods and 
data sets used in (3 determinations from peculiar veloc- 
ity studies. The IRAS result is tighter than our result 
because <7 8 ./ is better constrained and more peculiar ve- 
locity data have been compared to IRAS predictions. 
A weighted average of IRAS and 2MASS results gives 

cr 8 (ft m /0.3)°- 6 = 0.85 ±0.05. 

Feldman et al. (2003) measured the mean relative 
peculiar velocity for pairs of galaxies, and in compar- 
ison with the correlation function, derived a value of 
er 8 = 1.13 ± 0.22 for fl m — 0.30^' g£, which is consistent 
with our result. Feldman et al. do not quote a result for 
the combination f2 06 cr 8 , which one would expect to have 
smaller errors, than the marginal errors on ct 8 and Q m 
separately. 

Unlike the v-v comparisons, values can also be ob- 
tained by density-density (5 — S) comparisons. Previ- 
ously, this method required the construction of the grav- 
itation potential from the radial velocity field (POTENT; 
Bertschinger & Dekel 1989), differentiation and compar- 
ison it to the galaxy density field derived from redshift 
survey data. Note that such comparisons are based on 
smoothing and differentiating sparse and noisy peculiar 
velocity data. Using this method, Sigad et al. (1998) 
found (3\ — 0.89 ± 0.12. A newer method based on 
an unbiased variant of the Weiner filter yields values 
that are more consistent with those of the v-v method: 
(3i = 0.57 ± 0.12 (Zaroubi et al. 2002). 

It is possible to use peculiar velocities to obtain esti- 
mates of (3 in a completely different way. Peculiar veloc- 
ities of galaxies distort the pattern of galaxy clustering 
in redshift space, making the redshift space power spec- 
trum anisotropic. One can use the distortions to measure 
the parameter (3, and with information about the density 
field obtain the cosmological constraint osSl^ 6 . We com- 
bine the measured ct 8 g for the 2-degree Field Galaxy red- 
shift Survey (2dFGRS) by Lahav et al. (2002) with the 
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TABLE 3 

A SUMMARY OF THE (3 RESULTS FROM THIS PAPER AND RECENT V-V DETERMINATIONS OBTAINED FROM THE LITERATURE 



Comparison Sets 



<7 8 



Reference 



SN la vs. 2MASS 
SBF vs. 2MASS 
SFI vs. 2MASS 
SN Ia+SBF+SFI vs. 



2MASS 



0.47 ± 0.06 
0.46 ±0.06 
0.55 ±0.05 
0.49 ±0.04 



0.86 ±0.15 
0.84 ±0.15 
1.02 ±0.15 
0.91 ±0.12 



This Study 
This Study 
This Study 
This Study 



Mark III vs. IRAS 1.2 Jy 0.50 ±0.10 0.82 ±0.16 

SN la vs. IRAS 1.2 Jy 0.40 ±0.15 0.65 ± 0.25 

SBF vs. IRAS 1.2 Jy 0.42±°;£° 0.70±±g;^g 

SFI vs. IRAS 1.2 Jy 0.60 ±0.10 0.99 ± 0.09 

Mark III vs. IRAS 1.2 Jy 0.50 ± 0.04 0.82 ± 0.08 

Mark III vs. PSCz 0.60 ±0.10 0.99 ±0.19 

ENEAR vs. PSCz 0.50 ±0.10 0.82 ±0.16 

SN la vs. PSCz 0.54 ±0.06 0.89 ±0.10 

SFI vs. PSCz 0.42 ±0.04 0.66 ± 0.08 

SEcat vs. PSCz 0.51 ± 0.06 0.84 ± 0.10 

Weighted IRAS Avg. 0.50 ± 0.02 0.83 ± 0.06 



Davis ct al. (1996) 
Riess et al. (1997) 
Blakeslee et al. (1999) 
da Costa ct al. (1998) 
Willick & Strauss (1998) 
Saunders ct al. (1999) 
Nusser et al. (2001) 
Radburn-Smith et al. (2004) 
Branchini ct al. (2001) 
Zaroubi et al. (2002) 



2MASS and IRAS Avg. 



0.85 ± 0.05 



2dFGRS (3 result of Hawkins ct al. (2003) to obtain the 
combination cr 8 (f7 m /0.3) ' 6 = 0.89±0.23. Similarly, Per- 
cival ct al. (2004) used a spherical harmonics method on 
the 2dFGRS to find cr 8 (Sl m /0.3) a6 = 0.95 ± 0.12. These 
are in good agreement with our result. 

The comparisons between peculiar velocity derived val- 
ues of cr 8 f2 m 6 can be extended to other cosmological 
data. For instance, Refregier (2003) compiled an average 
o"8 = 0.83±0.04 value based on many weak lensing survey 
results for the assumed cosmology defined by S! m = 0.3, 
Q A = 0.7, T = 0.21. Contaldi et al. (2003) used Cos- 
mic Microwave Background (CMB) plus weak lensing re- 
sults to constrain the combination at cr 8 (f2 m /0.3) ' 5 = 
0.88 ±0.04. Tegmark et al. (2004) used recent CMB and 
Sloan Digital Sky Survey (SDSS) power spectra to find 
as (rim/0.3) 0,6 = 0.95 ± 0.13 under the assumption of a 
flat, n = 1 universe. Similarly, Seljak et al. (2005), com- 
bined WMAP with SDSS clustering, bias and Lya forest 
data to obtain (T8(f2 m /0.3) ' 6 = 0.86 ± 0.05 (assuming a 
flat universe with free n) . Using CMB and 2dFGRS data, 
Sanchez et al. (2005) find cx 8 (f2 m /0.3) a6 = 0.71 ± 0.10 
assuming (flat, n = 1). Sanchez et al. do not quote errors 
for the combination n 6 as, so above we have calculated 
the error assuming that the uncertainty in Sl m and as are 
independent. Since their Figure 1 indicates that these 
parameters are not independent, the above uncertainty 
will be an overestimate. Although the conflict is not sig- 
nificant, it is interesting that the Sanchez et al. result 
is lower than the CMB plus SDSS values, the 2dFGRS 
redshift-space distortion results, as well as our result. 

Overall, our peculiar velocity results are in good agree- 
ment with a broad variety of recent 08f2 m 6 measurements 
from the literature. 

5.3. Implications for Large-Scale Flows 

Although the focus of this paper has been on the grav- 
ity and peculiar velocity field within i? max , there are im- 
portant gravitational sources beyond this limit. While 
the existence of these sources do not affect our (3k re- 
sults (as we have preformed the analysis in the LG frame 
where external contributions are negligible), we can turn 



the problem around and place constraints on the large 
scale flow using our value of [3 k- 

Specifically, for the best fitting (3 K = 0.49 ± 0.04, we 
calculate a LG velocity v(0) = 403 ± 27 km s^ 1 in the 
direction I = 257°, b = 37°. The uncertainty on this 
linear theory predict ion is underestimated, because we 
expect that part of the LG's motion is not well described 
by linear theory. Since the LG is a galaxy group, and 
the relative infall of the Milky Way and Andromeda has 
already been accounted for following Yahil et al. (1977), 
it should be less affected by nonlinearities than are indi- 
vidual galaxies. It seems reasonable to adopt a thermal 
velocity dispersion of 100 km s , making v(o) uncertain 
at the 25% level. Thus for our best fit (3, the local volume 
does not fully account for the LG's motion with respect 
to the CMB, and there is a residual velocity dipole of 
U = 271 ±104 km s" 1 in the direction^ 300°, b = 15°. 

This residual dipole is slightly lower than, but consis- 
tent with, the result of Hudson (1994b), who compared 
the predictions of an optically-selected galaxy density 
field to D n — a and Tully-Fishcr peculiar velocities and 
found (3 = 0.50 ± 0.06 and a residual velocity, from be- 
yond SO/i^Mpc, of 405 ± 45 km s" 1 toward I = 292°, 
6 = 7°. The residual LG dipole found here is in good 
agreement with a recent estimate of the bulk flow of 
U = 225 km s" 1 toward I = 300°, b = 10° found by 
Hudson et al. (2004), for a sample of peculiar velocities 
at depths greater than 60/i _1 Mpc. 

Thus, our derived (3k suggests that there are significant 
contributions (<~ 40 ± 15%) to the LG's motion arising 
from sources at r > 65/i _1 Mpc. 

6. CONCLUSIONS 

In this paper we have used the 2MASS catalog and 
public rcdshift data to reconstruct the local density field. 
Under the assumption of GI and we used linear theory 
to derive the peculiar velocity fields which we compared 
to measured peculiar velocity data. The different data 
sets all yield results from the VELMOD method that 
arc in very good agreement, with a best fit (3 — 0.49 ± 
0.04. Calculation of the r.m.s. density fluctuations of 
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K-selected galaxies allowed us to generalize the result as 
a 8 (O m /0.3) - 6 = 0.91 ± 0.12. Our result is consistent 
with previous determinations of redshift space distortion 
methods, CMB and weak lensing results. The low value 
of (3 suggests that significant contributions to the LG's 
motion arise from beyond 65/i _1 Mpc. 

We would like to stress the power of the method 
that we have outlined here. Peculiar velocity compar- 
isons provide direct and independent measures of agfl^f, 
which are consistent with results from a wide range of 
techniques. Peculiar velocities probe scales near the fidu- 
cial 8/i _1 Mpc scale and allows us to make these pre- 
dictions without any assumptions regarding the back- 
ground cosmology. All that is required is a source of red- 
shift data, and a sample of galaxies which have redshift- 
independent distance information. This illustrates the 
power of peculiar velocity studies to probe the underly- 
ing cosmology. 

In future peculiar velocity work, it may be possible to 
break the degeneracy between f2°; 6 and b by applying 
a more sophisticated biasing model, such as a nonlinear 
bias based on the halo model (e.g. Marinoni & Hudson 
2002). Improving the predictions in the mildly nonlinear 



regime (e.g. Nusser & Branchini 2000) will be important 
for this latter goal. 

Furthermore, a deeper and more complete 2MASS data 
set will reduce the errors associated with sparseness and 
enable a better reconstruction of the local density field. 
Similar improvement will come with more reliable dis- 
tance estimates and larger peculiar velocity surveys. For 
instance, Wang et al. (2005) are currently working on 
compiling a sample of SN la that have distance errors, 
as low as 3 — 4%, although the sample is still quite small 
at present. As the data and techniques continue to im- 
prove, the future of peculiar velocity studies will continue 
to be a promising area for understanding the mass dis- 
tribution in the local universe. 
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